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Abstract 

The  formation  of  damage,  which  results  from  the  large  volume  expansion  of  the  active  sites  during  electrochemical  cycling,  in  rechargeable 
Li-batteries,  is  modelled  from  a  fracture  mechanics  viewpoint  to  facilitate  the  selection  of  the  most  effective  electrode  materials  and  configu¬ 
rations.  The  present  study  is  a  first  step  towards  examining  stable  cracking  in  such  high-energy  storage  devices,  by  considering  three  different 
configurations  at  the  nanoscale,  which  are  currently  at  an  experimental  stage.  As  a  result,  stability  diagrams  concerning  crack  growth  are 
constructed  and  compared  for  the  following  cases:  (a)  the  electrodes  are  thin  films,  (b)  the  Li-insertion  sites  in  the  anode  are  nanofibre-like 
inclusions,  (c)  the  active  sites  in  both  electrodes  are  spherical. 
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1.  Introduction 

Due  to  the  small  Li-intercalation  of  carbon,  which  is  used 
as  a  base  material  for  negative  electrodes  in  rechargeable  Li- 
batteries,  extensive  research  is  being  performed,  for  over  two 
decades  in  order  to  find  alternative  anode  materials.  This  re¬ 
search  has  suggested  that  some  of  the  best  candidates  are  Sn 
and  Si,  due  to  their  high  capacity  (990  and  4000  mAh  g_1). 
These  materials,  however,  have  not  been  used  commercially 
because  of  their  large  volume  expansion  upon  Li-insertion 
(during  the  charge/discharge  cycle),  which  results  in  crum¬ 
bling  and  severe  cracking  of  the  electrode  after  continuous 
electrochemical  cycling  [15].  It  is  therefore  anticipated  that 
modelling  fracture  in  these  high-energy  storage  devices  will 
provide  further  insight  towards  the  behavior  of  these  mate¬ 
rial  systems,  and  allow  for  the  selection  of  other  appropriate 
materials. 

A  first  attempt  to  model  the  mechanical  response  of  Li- 
batteries,  from  a  purely  elastic  point  of  view,  during  an  elec¬ 
trochemical  cycle  was  done  by  Aifantis  and  Hackney  [1].  The 
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stress-induced  inside  the  electrodes  was  modelled  by  treating 
them  as  thin  films  that  consisted  of  disc-shaped  Li-ion  active 
sites  embedded  in  an  inert  glass  or  ceramic  matrix  (Fig.  1). 
This  study  is  a  first  step  towards  modelling  the  most  important 
effect  that  this  cycling  deformation  process  has:  cracking.  As 
was  shown  in  experimental  evidence  provided  by  Aifantis  and 
Hackney  [1]  (Fig.  2),  continuous  electrochemical  cycling  on 
a  single  crystal  of  LfiVh^CL  resulted  in  multiple  fractures  on 
its  surface,  producing  nanocrystals  (Fig.  2).  It  is  believed  that 
continuous  cycling  would  have  the  same  effect  on  LiM^CL 
if  it  were  the  active  site  of  the  cathode;  the  individual  nanopar¬ 
ticles,  however,  produced  by  fracture,  would  no  longer  be  in 
electrical  contact  with  the  remainder  of  the  electrode.  As  a 
result  the  material  that  fractures  loses  electrical  contact  with 
the  electrode  and  becomes  unable  to  respond  to  the  applied 
voltage  required  for  charging  or  controlling  the  discharge  of 
the  battery;  hence  the  electrode  becomes  unsuitable  for  fur¬ 
ther  use  [1,2]. 

Moreover,  the  fracture  of  individual  particles  increases 
the  surface  area  available  to  chemical  attack  by  the  corro¬ 
sive  agents  of  the  battery  (HF  and  residual  H2O)  that  are 
believed  to  attack  the  surfaces  of  the  active  material  [1].  In 
fact,  it  is  reasonable  to  expect  that  the  chemical  instability 
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Fig.  1.  Idealized  geometry  of  the  electrode:  Li-insertion  particles  (shaded) 
embedded  in  a  glass  (blank)  matrix.  A  unit  cell  is  defined  by  a  circle  of  radius 
b  surrounding  a  circular  particle  of  radius  a. 

of  the  active  particle  surfaces  may  interact  with  the  mechan¬ 
ical  stress  to  enhance  the  structural  instability  of  the  material 
component  under  consideration  (stress  corrosion  cracking). 
Many  battery  developers  have  taken  the  approach  of  using 
large  particle  sizes  in  order  to  reduce  the  surface  area  avail¬ 
able  to  chemical  attack.  However,  this  will  accentuate  the 
problem  associated  with  stress  concentrations  at  the  surface 
of  the  particles.  An  alternative  method  that  can  be  used  is 
to  develop  a  composite  material  in  which  nanoscale  elec- 
trochemically  active  material  is  encased  in  large  particles 
of  a  chemically  inert  matrix  [3].  This  method  not  only  re¬ 
duces  the  surface  area  of  the  active  material  available  to 
chemical  attack,  but  also  minimizes  the  gradients  in  con¬ 
centration  responsible  for  fracture,  since  the  active  material 
component  of  the  composite  has  a  nanometre  length  scale. 
Therefore,  throughout  the  present  work,  both  electrodes  are 
taken  to  consist  of  Li-active  sites  (with  a  nanometre  diam¬ 
eter)  embedded  periodically  in  an  inert,  with  respect  to  Li, 
matrix. 

Consideration  of  the  above  electrode  configuration  along 
with  fabrication  methods  that  are  in  the  experimental  stage 
allows  the  mechanical  modelling  for  the  following  three 
cases:  (a)  the  active  sites  are  disc-shaped  platelets  and 
hence  the  electrodes  take  the  form  of  thin  films;  (b)  the  ac¬ 
tive  sites  of  the  anode  are  fibre-like  inclusions  (long  cylin¬ 
ders);  (c)  the  active  sites  in  both  the  anode  and  cathode  are 
spherical. 


Fig.  2.  Schematical  representation  of  fracture  zone  in  a  single  crystal 
LiMr^CU  due  to  nanocrack  formation  (damaged  surface  layer). 


Fig.  3.  Radial  cracking  configuration  assumed  for  the  unit  cell. 


2.  Modelling  the  electrode  damage  zone 

Since  experiments  [1]  have  shown  that  significant  dam¬ 
age  results  from  continuous  electrochemical  cycling  it  can  be 
assumed  that  the  fracture  layer  surrounding  each  active  site 
(Fig.  2)  has  undergone  severe  fracture,  such  that  it  can  support 
only  radial  stresses.  The  damage  zone  can  thus,  be  modelled 
by  a  sufficient  number  of  radial  cracks  that  initiate  at  the  ac¬ 
tive  particle/matrix  interface,  as  shown  in  Fig.  3.  It  should  be 
noted  that  a  and  b  denote  the  radii  of  the  active  Li-insertion 
site  and  surrounding  matrix,  respectively.  The  pressure  the 
Li-ions  induce  into  the  active  site  is  taken  to  be  uniform,  and 
is  set  equal  to  p ,  while  the  pressure  at  the  outer  boundary 
of  the  unit  cell  is  taken  to  be  q.  It  should  be  noted  that  this 
configuration  is  valid  upon  maximum  Li-insertion,  therefore, 
A  denotes  the  maximum  expansion  of  the  active  sites  in  the 
absence  of  the  matrix.  Finally,  the  radii  of  the  radial  cracks 
that  form  as  a  result  of  the  continuous  charge/discharge  cycle 
are  assumed  to  be  equal  and  their  length  is  p  —  a.  Since  the 
damage  region  defined  by  their  length  p  —  a  may  be  viewed 
as  supporting  only  radial  stresses  (crr),  all  the  other  stress 
components  vanish  ( a q  =  a yq  =  0)  inside  this  zone.  This  ap¬ 
proach  is,  in  fact,  similar  to  that  adopted  by  Dempsey  et  al. 
[4]  in  addressing  a  similar  problem  in  ice  mechanics. 


3.  Thin  film  electrodes 

In  hope  of  increasing  process  efficiency  and  optimizing 
storage  capacity  (i.e.  Li-intercalation  in  the  negative  elec¬ 
trode)  significant  efforts  have  been  made  to  prepare  anodes 
that  are  comprised  from  thin  layers.  With  this  fabrication 
method  it  is  believed  that  it  will  be  possible  to  exceed  the 
electrochemical  capacity  of  bulk  anode  materials  by  increas¬ 
ing  the  amount  of  Li-insertion  sites  in  the  thin  film.  This  has 
been  successfully  performed  for  anodes  comprising  of  tin 


K.E.  Aifantis,  J.P.  Dempsey  /  Journal  of  Power  Sources  143  (2005)  203-211 


205 


composite  oxide  (TCO)  [5]  (it  should  be  noted  that  the  active 
site  is  the  Sn,  while  the  oxide  acts  as  the  inert  matrix).  An 
additional  feature  of  this  configuration  is  that  it  contributes 
to  the  miniaturization  of  such  devices  since  the  thickness  of 
these  films  can  be  1  nm. 

The  main  advantage,  however,  from  the  mechanics  point 
of  view,  is  that  this  thin  film  configuration  makes  it  feasible  to 
model  the  electrodes  assuming  a  state  of  plane  stress;  thus,  the 
Li-active  sites  are  treated  as  disc-shaped  platelets,  embedded 
in  the  inert  ceramic  or  glass  matrix.  Since  the  same  config¬ 
uration  (and,  hence,  electromechanical  mechanism)  may  be 
assumed  for  both  the  anode  and  cathode,  the  subsequent  anal¬ 
ysis  is  valid  for  both  electrodes  with  the  only  difference  being 
the  Li-insertion  site;  for  example,  LiM^CL  may  be  used  in 
the  cathode,  and  Sn  in  the  anode;  the  matrix  will  be  taken  to 
be  soda  glass,  according  to  [5]. 

3.1.  General  displacement  and  stress  expressions  for 
elasticity  in  polar  coordinates 

The  stress  and  displacement  expressions  in  the  undam¬ 
aged  active  site  and  matrix  can  be  modelled  by  using  isotropic 
linear  elasticity  in  polar  coordinates.  The  unit  cell  under  con¬ 
sideration  suggests  an  axially  symmetric  problem  (<jr6i  =  0); 
thus,  the  non-zero  stress  components  and  displacements  are 
found  by  considering  the  equilibrium  relation 

dar  1 

-L  + -(ar  -  cre)  =  0,  (1) 

or  r 

which  in  combination  with  the  Airy  stress  function,  Hooke’s 
law  and  strain  compatibility  [6]  gives 

A 

oy  =  —y  +  B{\  +  2 log  r)  +  2 C, 

yZ, 

<7q  =  — y  +  B(3  +  2  In  r)  +  2 C,  (2) 

yZ 


1  A(1  +  v) 

ur  =  — - 1-  2(1  —  v)Br  log  r  —  B(  1  +  v ) 

E  [_  r 

+  2C(1  —  v)r] ,  ue  =  -pyrO,  (3) 

E 

where  (r,  0)  denote  polar  coordinates,  ( E ,  v )  denote  the 
Young’s  modulus  and  Poisson’s  ratio,  and  the  constants  (A, 
B ,  C)  are  to  be  determined  from  appropriate  boundary  con¬ 
ditions. 

Since  the  displacement  uq  should  be  single- valued,  it  fol¬ 
lows  that  the  constant  B  vanishes  and,  thus,  the  appropriate 
expressions  for  the  relevant  stress  and  displacement  compo¬ 
nents  simplify  to 


A 


A 


ay  —  —r  +  2C,  <jq  —  — y  +  2C, 

yZ  yZ 


1 

Ur  =  ~ 


A(  1  +  v) 


+  2C(1  —  v)r 


(4) 

(5) 


As  was  mentioned  previously,  the  stress  inside  the  active 
site  is  uniform,  hence  As  is  zero  [this  also  follows  from  the 
mathematics,  since  A/(0)2  is  undefined]  and  Eqs.  (4)  and  (5) 
reduce  to 

oy  =  oq  =  2CS; 

2CS(1  -  vs)r 

Ur  =  - 

Es 

inside  the  Sn. 

3.2.  Stress  expressions  in  active  site  and  boundary 
conditions 

According  to  the  discussion,  pertaining  to  Fig.  3,  the 
stresses  at  the  interfaces  at  r  -  a  and  r-b  are  given  by 


(6) 

(7) 


crr(a)  =  -p,  crr(b)  =  -q.  (8) 

Combination  now  of  Eq.  (6)  and  the  first  expression  in  Eq. 
(8)  gives 

2CS  =  -p.  (9) 


Therefore,  the  stress  and  displacement  expressions  inside 
the  active  site  are  given  by 


<7r  =  CTQ  =  -p\ 


Ur 


-p(l  -  vs)r 

Es 


(10) 


Since  the  state  under  maximum  Li-insertion  is  being  ex¬ 
amined,  the  initial  stress-free  configuration  of  the  active  site 
is  taken  to  be  that  which  it  would  assume  upon  reaching 
its  maximum  expansion  in  the  absence  of  any  constraints 
imposed  by  the  surrounding  matrix.  Since  the  configuration 
is  radially  symmetric,  this  maximum  stress-free  volumetric 
expansion  of  the  Sn  corresponds  to  an  increase  of  its  ra¬ 
dius,  denoted  by  A;  therefore,  the  initial  radius  of  this  par¬ 
ticle  is  taken  to  be  r  =  a  +  A.  The  surrounding  matrix,  which 
is  present  under  the  given  confined  configuration,  opposes 
the  aforementioned  free  expansion  by  pushing  back  the  ac¬ 
tive  site  by  a  distance  8;  therefore,  the  final  radius  of  the 
Sn  is  r  =  a  +  A  —  8,  and  the  total  displacement  of  the  outer 
surface  of  the  Sn  particle  at  r  =  a  +  A  is  ur  =  —8.  Combina¬ 
tion  of  this  condition  with  the  second  expression  in  Eq.  (10) 
gives 


ur(a  +  A)  =  —  8  =  — 


vs)(a  +  A) 

Es 


(11) 


Now  the  displacement  at  the  active  site/matrix  interface 
needs  to  be  formulated.  The  initial  stress-free  configura¬ 
tion  of  the  matrix  annulus  is  taken  to  be  that  at  which 
no  Li-ions  have  diffused  inside  the  active  sites,  so  the  ini¬ 
tial  inner  radius  of  the  glass  is  r  =  a.  As  the  inner  sur¬ 
face  of  the  glass  annulus  and  the  outer  surface  of  the  Sn 
particle  are  in  contact,  it  follows  from  the  above  para¬ 
graph  that  once  the  Sn  reaches  its  maximum  expansion  A, 
the  matrix  pushes  back  a  distance  8.  Thus,  the  final  inner 
radius  of  the  glass  is  r  =  a  +  A  —  8,  and  the  displacement 
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condition  for  the  glass  at  the  Sn/glass  interface  ( r  =  a )  is 
ur  =  A  —  8.  Combining  this  boundary  condition  with  Eq.  (11) 
gives 


ur(a)  —  A  —  8  —  A  — 


vs)(a  +  A) 
~E~s 


(12) 


3.3.  Stress/displacement  expressions  in  the  matrix 


case,  it  is  modified  to  plane  stress  ([8],  p.  103),  as 

(b2/r2)  +  ( (1  -  ug)/(l  +  ug)) 


u+(r)  = 


2pi, 


q 


p* 


(p2/r2)  +  { (1  -  Wg)/(1  +  Wg)) 
1  — (p2/b2) 


(b2/ p2)  —  1 


for  p  <  r  <  b, 


(18) 


The  corresponding  stress  and  displacement  expressions 
within  the  damage  region  are  determined  by  following  an 
analysis  similar  to  that  done  by  Dempsey  et  al.  [4] .  This  dam¬ 
age  zone  is  modelled  as  a  system  of  radial  cracks,  such  that 
a  uniaxial  state  of  stress  is  assumed  to  exist  within  it,  i.e. 
cf0  =  crro  =  0.  Thus,  the  equilibrium  relation  given  in  Eq.  (1) 
reduces  to 

d<7r  or  k 

— - 1 - =  0  =4  crr(r)  =  -,  for  a  <  r  <  p,  (13) 

dr  r  r 

where  k  is  an  integration  constant.  Given  that  ar(a)  =  —p,  it 
follows  that  k-  —  pa.  Furthermore,  Hooke’s  law  gives 


(7  r  —  E  i 


d  ur 
dr 


d  ur 
dr 


pa 


E0r 


urir ) 


Pa  !  ,  .  ,  * 

ln(r)  +  u  , 

Eo 


for  a  <  r  <  p, 


(14) 


where  an  expression  for  the  constant  u  can  be  found  by 
setting  the  displacement,  ur(p ),  right  in  front  of  the  crack 
tip  (i.e.  just  inside  the  uncracked  region)  equal  to  a  constant 
u+(p )  and  then  substituting  back  into  Eq.  (14) 

pa 

ur(p)  =  u+(p)  =>  u+(p)  =  —  —  ln(p)  +  u  =E  u 

Eg 

pa 

=  —  ln(p)  +  w+(p),  (15) 


then  insertion  of  Eq.  (15)  in  Eq.  (14)  concludes  that  inside 
the  cracked  region  the  radial  displacement  is  given  by 


ur(r)  = 


for  a  <  r  <  p. 


(16) 


The  uncracked  region  lies  between  the  crack  tip  bound¬ 
ary  (r=  p)  and  the  glass/glass  interface  ( r-b ),  and  it  can  be 
treated  as  a  hollow  disc  that  is  subject  to  an  internal  pressure 
p*  and  an  external  pressure  q  (exerted  from  the  neighbouring 
unit  cell).  The  internal  pressure  is  equal  to  that  present  at  the 
interface  with  the  fractured  region  (i.e.  at  r-  p)  and  it  can  be 
found  by  direct  substitution  in  Eq.  (13),  i.e. 


pa 

P *  =  -(Jrip)  =  — • 

P 


(17) 


The  displacement  solution  for  a  plane  strain  configuration 
(i.e.  a  hollow  cylinder)  that  is  subjected  to  an  internal  and 
external  pressure,  is  given  by  Westergaard  [7] ;  for  the  present 


where  /xg  =  Eg/ [2(1  +  ng)]  is  the  shear  modulus  of  the  un¬ 
cracked  glass  matrix  and  u+(r)  denotes  the  corresponding 
displacement  in  this  region. 

Thus,  by  setting  r=  p  in  Eq.  (18),  the  following  expression 
for  u+(p )  is  obtained 


u+(p)  = 


p[b\p*  +  P*V g  -  2 q)  -  p2p*(v g  -  1)] 
2p(vn  +  1  )(b2  -  p2) 


(19) 


Now,  insertion  of  Eq.  (19)  in  Eq.  (16)  gives  a  second  ex¬ 
pression  for  the  displacement  at  r  =  a 


ur{a) 


pa 

E 


b2  —  Cbp 

4-  2 _ - 

^  b2-p2 


(20) 


where  C  =  qb/pa.  The  displacement  at  r-b  can  also  be  de¬ 
duced  from  Eq.  (18)  as 


ur(b)  ~ 


pa 


bp  —  Cb 2 

2  b2  _  p2  +  C(1  +  Ug) 


(21) 


3.4.  Stability  index  formulation 


The  question  of  cracking  in  rechargeable  Li-batteries  may 
now  be  addressed  by  considering  the  hoop  stress  (ere)  just 
outside  the  damage  zone.  This  is  so  because,  for  the  present 
configuration,  gq  may  be  thought  of  as  being  the  opening 
tensile  stress  responsible  for  crack  stability  and  growth.  Using 
familiar  concepts  from  linear  elastic  fracture  mechanics  for 
radial  symmetry  [9]  the  energy  release  rate,  G,  can  be  defined 
as  a  function  of  crack  length  for  this  configuration  as 


G(p)  =  f3-a2(p+), 
Eg  n 


(22) 


where  n  is  the  number  of  radial  cracks,  and  oo{p+)  can  be 
deduced  from  [4]  as 


,  pa  11  +  ip/b)[ip/b)  -  2 C] 

<70{p  )  =  -  - ■= -  •  . 

b  l  (p/b)[l  -  (p/b)2]  J 

Now,  following  [4]  the  stability  index  is  defined  as 
b  d  G 

K  = - . 

G  dp 


(23) 


(24) 


It  can  be  seen  that  the  energy  release  rate  depends  on 
the  material  parameters  (E,  v),  the  geometric  parameters  ( a , 
b ),  the  number  of  cracks  n ,  as  well  as  on  the  internal  and 
external  pressures  p  and  q.  The  internal  pressure  can  be  found 
by  equating  Eqs.  (20)  and  (12)  and  solving  for  p,  while  the 
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external  pressure  q  need  not  be  defined  explicitly  as  long  as 
C  is.  Below  are  solutions  for  three  different  outer  boundary 
conditions. 


3.5.  Three  different  outer  boundary  conditions 


3.5.1.  Case  1:  clamped  outer  boundary 

The  first  case,  which  is  considered,  corresponds  to  the 
common  “manufacturing  consistent”  configuration  accord¬ 
ing  to  which  the  battery  system  is  tightly  constrained  by 
the  outer  casing  and,  therefore,  the  displacement  at  the 
glass/glass  interface  is  zero.  Thus,  ur(b)  =  0  and  in  view  of 
Eq.  (21),  one  can  solve  for  C  to  obtain 


2  bp 

(1  -  vg)b2  +  (1  +  ug )p2' 


(25) 


Then  inserting  Eq.  (25)  in  Eq.  (20),  and  setting  the  result¬ 
ing  expression  equal  to  Eq.  (12),  the  internal  pressure  can  be 
explicitly  calculated  as  a  function  of  the  geometric  and  ma¬ 
terial  parameters.  Finally,  combination  of  Eq.  (22)  with  Eq. 
(23)  yields  the  energy  release  rate  for  this  case,  G\ ,  as 


Gi(p)  = 


where 


p\  =  A 


2  2 
np^a 

nEgp 


(1  -  vg)b2  -  (1  +  vg)p2 
(1  -  vOb2  +  (1  +  Vo )p2 


a 


2(1  -  vg  )b‘ 


(1  -  Va)b2  +  (1  +  Va)p2 


-  (1  -  Vo) 


+ 


{a  +  A) 

f 


-l 


(26) 


(27) 


It  should  be  noted  that  in  Eq.  (27),  as  well  as  in  the  ex¬ 
pressions  to  follow,  we  define  T  =  Es/(  1  —  vs). 


3.5.2.  Case  2:  traction  free  outer  boundary 

The  second  outer  boundary  condition  was  formulated  by 
considering  the  “natural”  condition,  according  to  which  the 
pressure  that  is  induced  on  the  glass  by  the  active  site  fades 
with  increasing  distance  and  hence  the  outer  pressure  (at  the 
glass/glass  interface)  is  zero.  Therefore,  q  =  0  (which  implies 
that  C  =  0),  and  thus,  Eqs.  (20)  and  (12)  can  be  used  as  before 
to  determine  the  internal  pressure,  for  this  case,  p2 .  Then  by 
using  Eqs.  (22)  and  (23)  the  energy  release  rate,  G2,  is  readily 
calculated  as 


G2(P)  = 


Ttp2a2  b2  +  px 


nEgp 


b2  —  p2 


(28) 


where 


P2  =  A 


fW(-)+2 

Eg  |_  \al 


bx 


b2  —  p2 


-  (1  -  Va) 


+ 


(1 0  +  A) 

f 


(29) 


3.5.3.  Case  3:  self-equilibrated  loading 

The  final  case  to  be  considered  is  that  of  “self-equilibrated 
loading”,  according  to  which  the  force  (2 j;qb)  that  is  exerted 
on  the  glass  annulus  by  the  surrounding  unit  cell  is  equal  to 
that  exerted  onto  it  by  the  Li-insertion  site  ( l^upa ).  Therefore, 
C  =  1 ,  and  solving  for  the  internal  pressure  as  before,  the 
energy  release  rate,  G3 ,  for  this  case  is  found  to  be 


Gffp)  = 


7tp2a2  fb  —  p 
nEgp  \b  +  p 


2 


(30) 


where 


(31) 


3.6.  Stable  crack  growth  in  negative  electrode 

Finally,  the  stability  index,  /c,  can  be  computed  for  each 
case,  via  Eq.  (24),  by  letting  the  radii  ( a  and  b)  and  the 
material  parameters  (Es,  Eg,  vs,  vg)  be  (100  and  lOOOnm) 
and  (41,  75GPa,  0.33,  0.25)  [10],  respectively,  while  the 
misfit  parameter  A  for  these  geometric  parameters  has  been 
found  to  be  145  nm  [1].  (It  should  be  noted  that  n  need  not 
be  quantified  since  k  is  independent  of  the  number  of  radial 
cracks  that  form.)  Now  a  stability  diagram  can  be  constructed 
by  plotting  k  with  respect  to  p  (Fig.  4).  Stable  crack  growth 
is  examined  with  the  energy  release  rate  G  being  exactly 
matched  by  the  resistance  to  fracture  R.  An  additional  nec¬ 
essary  condition  is  that  the  rate  of  increase,  with  change  in 
crack  length,  of  G  be  less  than  or  equal  to  the  rate  of  increase 
of  R.  Quasibrittle  materials  require  the  stability  index  k  to 
be  negative  for  stable  crack  growth;  it  follows  that  the  more 
negative  the  value  of  /c,  the  more  stable  the  crack  growth 
becomes,  since  the  difference  between  final  and  initial  G 
increases.  In  this  paper,  it  will  be  assumed  that  stable  crack 
growth  (slow  crack  propagation)  is  feasible  only  if  k  is  nega¬ 
tive.  As  the  stability  index  increases,  crack  growth  becomes 


Fig.  4.  Stable  crack  growth  in  negative  electrode  for  thin  film  (plane  stress) 
configuration. 
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more  and  more  unstable.  Once  the  p- axis  is  crossed,  unstable 
crack  growth  ensues  (catastrophic  failure  is  approached  very 
rapidly),  hence  the  p-intercept  indicates  the  critical  distance 
to  which  a  crack  can  propagate  stably  (the  stable  critical 
crack  length  is  found  by  subtracting  the  p-intercept  from 
the  inclusion  radius  a).  It  should  be  noted  that  in  the  case 
that  an  asymptote  is  featured  in  the  stability  diagram  (as  in 
Fig.  4),  the  crack  will  close  shut  well  before  the  asymptote, 
p-axis  intercept,  is  reached.  In  this  case,  the  energy  release 
rate  is  decreasing  with  crack  growth;  at  the  asymptote  G  =  0. 
In  particular,  it  can  be  seen  from  Fig.  4  that  for  Case  2,  the 
corresponding  critical  stable  crack  length  is  rather  small, 
since  the  p-axis  is  crossed  at  p  =  450  nm,  while  for  Case  1 
(C  >  1)  an  asymptote  is  attained  at  p  =  755  nm,  which  implies 
that  the  crack  will  close  shut  before  this  distance  is  reached. 
Finally,  for  Case  3,  the  stability  diagram  suggests  that  stable 
crack  growth  takes  place  until  the  outer  radius  b. 

It  can,  thus  be  concluded  that  the  “manufacturing  consis¬ 
tent”  case  is  more  desirable  since  the  cracks  are  pinched  shut 
once  their  critical  stable  length  is  approached. 


4.  Fibre-like  active  sites 


Another  fabrication  method  that  is  being  employed  to  in¬ 
crease  the  Li-intercalation  in  the  anode  is  that  according  to 
which  the  active  sites  take  the  form  of  fibres,  whose  dimen¬ 
sions  are  in  the  sub-micron  scale.  Experimental  observations 
showed  that  this  method  proved  to  be  succesful  in  the  case 
where  the  active  material  was  silicon  [11]  or  carbon  [12], 
since  this  allowed  a  higher  metal  content  to  enter  the  Li- 
insertion  site.  Other  anode  materials  which  are  still  in  the  ex¬ 
perimental  stage  are  Li-capacity  templated  anodes  consisting 
of  about  100  nm  Sn02  nanofibres  [1,13]. 

To  analyse  the  behaviour  of  the  aforementioned  anodes, 
their  configuration  is  taken  to  be  that  of  long  cylindrical  ac¬ 
tive  site  inclusions  embedded  in  an  inert  matrix,  and  thus, 
conditions  of  plane  strain  can  be  used.  The  plane  stress  anal¬ 
ysis  in  Section  3  can  be  modified  for  plane  strain  by  making 
the  following  substitutions  in  Eqs.  (1)— (3 1)  ([8],  p.  103), 


E 

v 


E'  = 


v  = 


1  — u2 

V 

l—v 


r 


rr  = 


(1  _  2us)(l  +  vs) 


(32) 


Fig.  5.  Stable  crack  growth  in  negative  electrode  with  nanohbre-like  active 
inclusions  (plane  strain). 


5.  Spherical  active  sites 


As  was  mentioned  in  the  introduction  the  third  case  that 
shall  be  examined  is  that  according  to  which  the  electrodes 
are  comprised  of  spherical  active  sites  embedded  in  a  matrix. 
Studies  concerning  this  electrode  configuration  are  being  per¬ 
formed  for  nanocomposite  Sn/C  [14],  Sn-Fe-C  [15]  as  well 
as  for  TCO  [1]  anodes.  In  all  these  systems  the  Sn  acts  as  the 
active  material. 

The  quantities  considered  are  the  same  as  in  those  in  Fig.  3, 
but  since  we  are  in  three  dimensions  now  the  stress  and  dis¬ 
placement  relations  change  accordingly  to 


2C  1  +  v 
&rr  —  y  T  2-  ——D\ 

1  —  2v 

C  1  +  v 

aee  —  Off  —  — t  T  2- - —  D\ 

1  —  2v 


ur  = 


2(1  +  v)  (  C 


2  r2 


+  Dr 


(33) 


while  all  other  stress  (crr0,  orQ,  o^q)  and  displacement  ( uq , 
u^)  components  vanish  due  to  symmetry. 

Again,  the  stress  inside  the  spherical  Sn  sites  does  not  vary 
with  position  and  hence  Cs  =  0.  Thus,  Eq.  (33)  simplifies  to 


Grr 


1  T 

=  &00  —  ^(jxp  =  2-  —  Ds , 

1  2vs 


U  y 


2(1  +  v)r 
E~s 


(34) 


Therefore,  similarly  as  before  a  stability  diagram  can  be  con¬ 
structed  for  this  case  as  shown  in  Fig.  5.  According  to  [13], 
the  active  sites  in  the  anode  are  taken  to  be  Sn,  and  the  matrix 
is  soda  glass. 

As  in  the  thin  film  configuration,  Case  2,  which  corre¬ 
sponds  to  “natural”  conditions  allows  stable  crack  growth  un¬ 
til  a  relatively  small  distance  (p  =  500  nm).  Under  the  “man¬ 
ufacturing  conditions”,  Case  1,  the  crack  will  close  once  the 
stable  critical  crack  distance  p  =  7 10  nm  is  approached,  while 
the  “self-equilibrated  loading”  case  (Case  3),  allows  for  sta¬ 
ble  cracking  until  the  outer  boundary  of  the  unit  cell. 


in  the  active  sites. 

Since  Eq.  (8)  still  holds,  the  expressions  in  Eq.  (34)  can 
be  re-written  in  terms  of  the  internal  pressure  p  as 

(1  —  2vs)  (1  —  2vs)rp 

arr  —  P  ^  —  p  Ur  — 

2(1  +  Vs)  Es 

(35) 

The  physical  reasoning  concerning  the  contrac¬ 
tion/expansion  of  the  Sn  and  surrounding  matrix  is  the 
same  as  in  Section  3,  but  in  view  of  the  new  stress  and 
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displacement  expressions  established  above,  Eq.  (11)  takes 
the  form 


u,(a  +  A)  =  —  8  =  — 


(1  -  2 vs)(r  +  a)p 


(36) 


while  the  displacement  at  the  Sn-glass  interface,  Eq.  (12),  is 
modified  to 


ur(a)  =  A  —  8  =  A  — 


(1  -  2 vs)(A  +a)p 


(37) 


As  was  already  stated  the  cracked  region  supports  only 
radial  stresses;  hence  the  force  equilibrium  relations,  under 
spherical  symmetry  reduce  to 


dov  2ar 

+  —  =0 


(38) 


dr  r 

solution  of  which  provides  the  stress  inside  the  damage  zone 


k 

cFrrif)  =  -y,  for  a  <  r  <  p, 


(39) 


where  k=  —pa2,  since  the  first  expression  in  Eq.  (8)  implies 
that  <7rr(a )  =  —p  =  k/r2 . 

Since  the  other  stress  components  inside  this  damage 
zone  vanish,  the  displacement  expression  can  be  determined 
through  the  relation  of  arr  with  ur  as 


^rrif)  —  E'g 1 


d  Ur 
dr 


d  ur 
dr 


pa 

HAr2 


,  x  P°  ,  * 

u,(r)  =  — - h  u  , 


Ear 


(40) 


where  the  constant  u  is  found  by  continuity  of  displacements, 
as  in  Section  3.  Thus,  ur(p)  is  set  equal  to  the  displacement 
right  in  front  of  the  crack  tip  u+(p)  to  give 


,  v  ,  v  Pa ,  * 

U  +  (p)  =  Ur(p)  =  — - h  u 


U*  =  u  +  (p)  — 


pax 


(41) 


Egp  Eg p 

It  follows  that  inserting  Eq.  (41)  in  Eq.  (40)  provides  the 
displacement  inside  the  cracked  region 

2  /l  1 


ur(r)  = 


pa 


- )  +  U+(P )» 

r  p 


for  a  <  r  <  p,  (42) 


while  the  displacement  inside  the  uncracked  glass  region  (i.e. 
for  a  hollow  sphere  subjected  to  an  internal  and  external  pres¬ 
sure)  is  given  by  Westergaard  [7]  as 


u+(r)  = 


r(  1  +  vg) 


p* 


bfjpPp  +  (1  —  2vg)(l  +  Dg) 

b2  /  p2  —  1 


-q 


p3/(2r3)  +  (1  -  2ug)(l  +  ve) 


1  —  p2/b 2 


p  <  r  <  b 


(43) 


where  p*  is  the  pressure  exerted  at  r=  p,  and  is  found  through 
Eq.  (39)  to  be 

.2 


P*  = 


pa 

P: 


(44) 


Now  the  displacement  right  in  front  of  the  crack  tip  can 
be  computed  by  direct  substitution  in  Eq.  (43) 


w+(/>)  = 


P 


2Eg{b2  -  p2) 
—q[3b2(l  -  vg)] J  . 


| /?*[(!  +  Vg)b2  +  2(1  -  2 Vg)p2] 


(45) 


Similarly  as  before  a  second  expression  for  the  displace¬ 
ment  at  the  Sn/glass  interface  ( r  =  a )  is  obtained  by  inserting 
Eq.  (45)  in  Eq.  (42) 


u(a)  = 


pa ' 


1 

a 


2(1  -  v8)  3(1  -  vg)(b3  -  Sbp 2) 


P 


2  p(b2  —  p2) 


(46) 


where  S  =  qb2/(pa2).  The  displacement  at  the  outer  glass 
boundary  ( r-b )  is  also  found  through  Eq.  (43) 


u(b)  = 


pa 

2bE< 


3(1  -  Vg)(b2p  -  Sb2) 
b2  —  p2 


+  <5(1  +  Vg) 


(47) 


Finally,  Eq.  (24)  for  the  stability  index  is  still  valid,  but 
in  spherical  coordinates  the  energy  release  rate  (according  to 
Dempsey  et  al.  [9])  is  defined  as 


2(1  -  Vg )pggg(p+) 

nE„ 


G(p)  = 


where, 


(48) 


m(p+)  = 


pa 

b2 


1  -  3  S(p/b)2  +  2 (p/bf 
2{p/b)\\  -  ( p/b )3) 


(49) 


5.7.  Three  different  outer  boundary  conditions 

5.1.1.  Case  1:  clamped  outer  boundary 

Since  u(b)  =  0  in  the  manufacturing  consistent  case,  Eq. 
(47)  is  set  equal  to  zero  and  S  is  found  to  be 


5  = 


3(1  -  vg)b2p 


2(1  -2iO£3  +  (l  +  tUp3' 


(50) 


Inserting  Eq.  (50)  in  Eq.  (46)  and  then  setting  the  result¬ 
ing  expression  equal  to  Eq.  (37),  gives  the  internal  pressure, 
which  can  in  turn  be  substituted  in  Eq.  (48)  to  give 


4^2 


G  i  = 


2(1  —  Vg)a  p 


nEgP' 


(1  -  2vg)b3  -  (1  +  ug)p3 
2(1  -2iv)63  +  (l  +iv)p3 


-i  2 


,  (51) 


where 


Pi  =  A 
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2  . 
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—  —  2(1  —  vg) 
pEg  la 

3(l-wg)(l-2  vg)b3 
2{\—2vg)bi+{\+Va)pi 
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(a  +  A) 

y^sph 


-1 


(52) 
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2x  10  7  4x  10  7  6x  10  7  8x  10  7  lx  10 

distance  p  (in  metres) 


Fig.  6.  Stable  crack  growth  in  negative  electrode  with  spherical  Sn  active 
sites  embedded  in  glass  matrix. 


2  x  1 0-7  4  x  1 0-7  6  x  1 0-7  8  x  1 0-7  1  x  10“6 

distance  p  (in  metres) 

Fig.  7.  Stable  crack  growth  in  negative  electrode  with  spherical  Sn  active 
sites  embedded  in  FeC  matrix. 


It  should  be  noted  that  in  Eq.  (52)  as  well  as  in  the  follow¬ 
ing  relations  Csph  =  Fg/(  1  —  2vg). 


5.1.2.  Case  2:  traction  free  outer  boundary 

For  the  traction  free  boundary  condition, 
<Jrr(b)  =  0  =>  q  =  0,  hence  S  =  0.  Inserting  this  value  of 
S  in  Eq.  (46)  and  then  setting  the  resulting  expression  equal 
to  Eq.  (37),  G2  is  found  as 


G2  = 


(1  —  Vg  )a4p2  [b3  +  2 p' 


-.2 


InEgp 3  \_b3  —  p3  \ 


where 


P2  =  A 


a 


pEi 


P  „„  N  ,  3(1  -  vg)b 3 
-  -  2(!  -  vg)  +  —3  - 


+ 


(a  +  A) 

y^sph 


-1 


(53) 


(54) 


5.1.3.  Case  3:  self- equilibrating  loading 

For  the  case  of  self-equilibrated  loading,  i.e.  pa2=qb 2, 
S  =  1  and  the  internal  pressure  is  solved  as  before  to  give 


G3  = 


(1  —  Vg  )a  p 


4„2  r 


(b  +  2p)(b  -  p) 


2nEgp3  |_  b2  +  bp  +  p2  J 


where 


P3  =  A 


ax 
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pEg 
(a  +  A) 

y^sph 


P 

a 


(1  -  Vg)  3(1  -  Vg )p‘ 


2  (b2  +  bp  +  p2) 


(55) 


(56) 


The  stability  diagrams  for  all  cases  considered  are  shown 
in  Fig.  6  (the  active  sites  are  again  taken  to  be  Sn,  according  to 
[14,15],  while  the  matrix  is  taken  to  be  soda  glass  for  consis¬ 
tency  with  the  previous  cases).  It  can  be  seen  that  Cases  1,  2 
and  3  allow  for  critical  stable  crack  growth  until  p-121  nm, 
560  nm  and  lOOOnm,  respectively. 


Since  Sn-Fe-C  composites  are  also  under  consideration, 
it  is  of  interest,  to  examine  crack  growth  when  the  matrix  is 
FeC,  as  [15]  suggests;  hence  Fig.  7  is  obtained  (p  =  680nm, 
600  nm  and  1000  nm  for  Cases  1,  2  and  3). 


6.  Conclusions 

The  present  study  is  a  first  attempt  to  predict  the  extent  of 
stable  crack  growth  in  Fi-ion  battery  electrodes  that  comprise 
of  active/inactive  nanocomposites.  In  particular,  stability  di¬ 
agrams  are  obtained  for  anodes  that  consist  of  Tin  Compos¬ 
ite  Oxide,  since  various  nano-configurations  (thin  film  elec¬ 
trodes,  fibre-like  active  sites  in  anode,  spherical  active  sites) 
still  in  the  experimental  stage  employ  this  material.  As  can  be 
seen  from  Figs.  4-7,  all  configurations  examined  herein  ex¬ 
hibit  similar  stability  behaviour  when  subjected  to  the  particu¬ 
lar  boundary  conditions  considered.  For  the  “manufacturing 
consistent”  case  it  is  desirable  to  have  small  critical  stable 
crack  lengths,  since  the  stability  diagrams  suggest  that  the 
smaller  the  critical  length,  the  smaller  the  distance  the  crack 
will  propagate  before  being  pinched  shut;  it  should  be  em¬ 
phasized  that  this  case  is  the  most  desirable  since  the  cracks 
cannot  propagate  beyond  a  certain  distance.  The  “natural” 
boundary  condition,  on  the  other  hand,  is  the  least  desirable 
since  it  is  inherently  unstable.  Finally,  no  instability  was  en¬ 
countered  for  the  “self-equilibrated  loading”  case;  the  crack 
could  propagate  stably  toward  the  outer  boundary  of  the  unit 
cell. 

It  follows  from  the  above  discussion  that  the  most  desir¬ 
able  configuration  of  all  those  considered  is  that  of  spheri¬ 
cal  Sn  sites  embedded  in  a  FeC  matrix,  having  not  only  the 
smallest  critical  stable  crack  length  for  the  “manufacturing 
consistent”  case,  but  also  the  largest  critical  stable  length  for 
the  “natural”  boundary  condition. 

A  direct  extension  of  the  present  work  is  to  apply  this 
analysis  to  other  material  systems,  such  as  those  which  use 
Si  as  the  active  site.  Moreover,  it  should  be  emphasized  that 
the  current  analysis  can  allow  for  the  determination  of  the 
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critical  crack  length  (pc)  at  which  fracture  of  the  electrode 
will  take  place  (by  substituting  the  critical  energy  release  rate 
of  the  matrix  in  Eqs.  (22)  or  (48)  and  solving  for  pc);  addi¬ 
tional  information  can  thus  be  obtained  as  to  which  materials, 
configurations  and  outer  boundary  conditions  are  more  ap¬ 
propriate.  It  may  be  possible,  for  example,  that  the  “natural” 
outer  boundary  condition  can  allow  for  greater  critical  crack 
lengths  (which  do  not  exceed  the  critical  stable  crack  length) 
than  the  “manufacturing  consistent”  condition. 

For  this,  however,  to  be  examined  the  critical  energy  re¬ 
lease  rate  for  soda  glass,  FeC  and  other  inert  ceramics  needs  to 
be  determined  through  experiments,  since  such  information 
is  lacking. 
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